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Abstract 



In a series of two recent papers, the frequency and size distribution dependence 
of extinction spectra for astronomical silicate and graphite grains was analyzed by us 
in the context of MRN type interstellar dust models. These grains were taken to be 
homogeneous spheres following the power law {a~^-^) size distribution which is very 
much in use. The analytic formulas were obtained for the graphite and silicate grains in 
wavelength range lOOOA - 22,500A and their utility was demonstrated. In this paper of 
the series, we present analytic formulas for the scattering and absorption spectrum of 
another important constituent of interstellar dust models, namely, the polycyclic aro- 
matic hydrocarbons (PAHs). Relative contribution of the PAHs to extinction vis a vis 
carbonaceous classical grains has been examined. 

Keywords: Interstellar grains. Extinction spectrum, PAH, Size distribution depen- 
dence. 



1 Introduction 



In a recent series of two papers, we presented analytic formulas for the frequency and size 
distribution dependence of the extinction spectrum of astronomical silicate and graphite 
grains in the wavelength range 1000 A - 22,500 A [1,2] (hereafter referred to as RSG-1 and 
RSG-2). These formulas were worked out keeping in mind the two component interstellar 
dust models of Mathis, Rumple and Nordsieck (MRN) [3] type, wherein one assumes the dust 
grains to be a collection of bare silicate and graphite homogeneous spheres, each of the two 
components obeying a power law size distribution with exponent —3.5, having a minimum 
radius (ao) and a maximum radius {am)- In our analysis, the admitted ranges of Oq and 
were 0.002 fim < aQ < 0.005 fim and 0.15 fim < < 0.25 fim for graphite grains and 
0.004 fim < ttQ < 0.006 fim and 0.2 fim < a„i < 0.4 fim for silicate grains. These ranges are 
thought to be quite appropriate and allow one to have the kind of flexibility needed in flxing 
the values of ag and a^- A comparison of extinction predicted by these formulas vis-a-vis 
exact Mie computations showed that these formulas give results which are in close proxim- 
ity to exact computations. Hence, these formulas have the potential to be used efficiently 
for assessment of the extinction contributions of the silicate and graphite components in 
investigations and building of MRN type models without resorting to large scale numerical 
computations. Moreover, by selecting proper ranges of frequency, the extinction behavior 
can be effectively approximated by yet simpler forms (e,g., linear, quadratic, etc.). That 
is, further simpliflcation can be made in an effort to set appropriate particle size ranges to 
start with. It may be mentioned here that analytic formulas fitting the extinction spectra of 
several stars have been obtained in the past by Cardelli et al. [4] and Fitzpatrick and Massa 
[5]. However, these parametrizations are mostly based on mathematical schemes. In con- 
trast, the formulas provided in RSG-1 and RSG-2, give some degree of physical insight into 
the problem and also address the problem in terms of silicate and graphite grains separately 
with a provision of allowable variations in oq as well as in - which are known to be quite 
adequate. 

In the MRN type models it is generally assumed that carbonaceous grains have ~ 
0.005 fim and their optical properties are more like graphite grains. However, in more recent 
models it has been recognized that the carbonaceous component also has a sub-component 
of "very small grains" whose optical properties are different from that of graphite. These 
grains are the so-called polycyclic aromatic hydrocarbons that are referred to as "astronom- 
ical" PAHs in the literature. 



1.1 Why PAH's 

In 1956 Piatt [6] proposed that very small grain or large molecules of radii around lOA to be 
present in the Interstellar Dust. Donn(1968)[7] proposed that polycyclic aromatic hydrocar- 
bon hke 'Piatt Particles' may produce UV extinction in the interstellar spectra. Importance 
and ubiquity of polycyclic aromatic hydrocarbons is quite evident in literature addressing 
the interstellar extinction. The inclusion of PAH's is based on the three decades of research 
by various groups. The IR emission features are very well attributed to PAH's by now and 
also allow to place constraints on size distribution of very small dust components. Ubiquity 
of PAH's in many galaxies has been observationally verified by ISO and Spitzer, thus mak- 
ing Polycychc Aromatic Hydrocarbons (PAH's) a very important constituents of Interstellar 
Dust(Draine and Li 2007,[8]). Draine and cowrokers (Li and Draine 2001b[9], 2002a[10]; 
Weingartener and Draine 2001a[ll]) proposed the Silicate- Graphite- PAH's model, of which 
PAH are the small-size end of the carbonaceous grains, is in excellent agreement with ob- 
servations. For more details please refer to the paper by Li and Greenberg[12]. On these 
lines, in this paper, we aim at extending the analyses of RSG 1 and RSG 2 to include the 
extinction spectra analysis of the PAHs in the size range 5 A < a < 50 A and spectral 
domain 1000 A to 22,500 A . The closed form analytic expressions for extinction by such 
grains should enable one to make more precise and elaborate model investigations and hence 
correct model building. The PAHs like grains can be further distinguished as being either 
in neutral or ionic state. Analytic formulas to evaluate and analyze the extinction spectrum 
have been obtained for both types of PAHs in this work. As the PAHs are very small in size, 
the scattering is expected to be unimportant [13]. Thus, the extinction can be effectively 
computed from the absorption only. However, for completeness, formulas have been obtained 
for absorption as well as scattering efficiencies. 

With this paper of the series, we aim to complete our objective of constructing an analytic 
framework by the use of which a direct approach to analysis of interstellar dust extinction 
spectra available for various galaxies viz. MW, LMC, SMC etc. can be done very ex- 
pediently. Whereas, it would be our endeavor to carry out the observed interstellar dust 
extinction spectra data analysis in context of MW, etc., in a forthcoming paper, neverthe- 
less, a comparison of the extinction produced by equal volumes of PAHs and the two varieties 
of graphite (perpendicular and parallel) separately is examined here to appreciate the im- 
portant role of PAHs in shaping the extinction spectra of the dust containing them. The 
mass density of PAHs and graphite is taken to be the same (?» 2.24 gm/cm^) here, as per 
the general convention so that equal volume also implies equal mass. 

The absorption and scattering efficiency data, covering a wide radius range 0.000355/xm < 
a < 0.01 jim and the wavelength range 10 A < A < 10^ A, for the "astronomical" PAHs 



in excellent tabular form can be found on the website of Draine [14]. The details in this 
regard are available in the appendix part of [9]. These results, to be parametrized in this 
paper, are referred to as the exact results. 

This paper has been organized as follows. Section 2 describes all the relevant functional 
forms of scattering and absorption efficiencies obtained for PAHs in the wavelength range 
lOOOA < A < 22, 500A for size range 5 A < a < 50 A. The accuracy of all the formulas 
obtained has been demonstrated by numerically computing the absorption and scattering 
efficiencies from the formulas and contrasting them with the corresponding exact results of 
Draine [14]. It may be mentioned that the parametrization of these so-called exact results in 
various wavelength domains have also been achieved by Li and Draine [9] . A typical numer- 
ical comparison of their formulas with exact absorption efficiencies has also been presented 
here in the UV and FUV regions. In section 3, the extinction spectra of various constituents 
of the carbonaceous component of the dust has been discussed and their relative contri- 
butions to the extinction spectrum have been analyzed. The convenience introduced by 
analytic formalism becomes very apparent in this part of the analysis. The dust model has 
been described contextually in this section. Finally, we conclude by summarizing and dis- 
cussing the salient results of this paper in section 4. In this work, grains are assumed to be 
spheres and, therefore, this work ignores issues related to polarization. 

2 Parametrization of PAHs absorption and scattering 
efficiencies 

The computed values of scattering and absorption cross sections of PAHs in the size range 
3.5 A < a < 100 A have been tabulated by Draine [14] for neutral as well as for ionized 
PAHs in the wavelength range lOA < A < lO^A . We have attempted to obtain sim- 
ple parametrization for corresponding absorption and scattering efficiencies for PAH sizes 5 
A < a < 50A for wavelength range 1000 A < A < 22,500 A. The PAHs in this size range 
have been taken to be the major contributors to the absorption (extinction) due to dust in 
Li and Draine [9] . Thus, it was decided to restrict the present investigation to the PAH size 
range 5A < a < 50 A. The admitted wavelength range is the same as it was in RSG-1 and 
RSG-2. 

For absorption efficiencies, the total wavelength range has been divided into five regions, 
corresponding to far-ultraviolet (1000 A < A 1800 A), ultraviolet (1800 A < A < 4000 A), 
visible (4000 A < A < 8000 A), infrared-I (8000 A < A < 12,500 A) and infrared-II (12,500 
A < A < 22,500 A). The corresponding formulas have been accomplished for the neutral as 
well as the ionic PAHs. Barring region 4 (infrared-I), in each of the other regions, it was 



possible to arrive at a single formula covering the entire size range 5 A < a < 50 A. In the 
region 4, separate formulas were required for the ranges 5 A < a < 10 A and 10 A < a < 
50 A. For the ionic PAHs however, separate formulas are required in regions 3, 4 and 5. 
In regions 1 (FUV) and 2 (UV), the corresponding absorption efficiency formulas for the 
neutral and the ionic particles are identical. For scattering efficiencies, we have divided the 
wavelength range of our interest in three spectral regions. Theses are far- ultraviolet (1000 
A < A < 1800 A), ultraviolet (1800 A < A < 4000 A) and visible- infrared (4000 A < A < 
22,500 A). 

A parametrization of absorption cross-section of PAHs has also been obtained by Li and 
Draine [9] starting from the far-ultraviolet to the far-infrared in the size domain 3.5A< a < 
lOOA. Their parametrization is characterized by a set of Drude profiles. Whereas for 
A <3300A the formulas require at most one Drude profile, as many as 12 Drude pro- 
files have been used to construct the formula for A > 3300 A. 



In this work, numerical comparisons have been generally presented for the values of ra- 
dius a = 10, 20, 35 and 50 A. But, it may be mentioned that, in addition, we have also 
verified all the results for the intermediate values of a = 5, 8, 15 and 45 A to ensure that 
all observations and statements made based on them are really speaking valid over the size 
range 5A < a <50A. 

2.1 Neutral PAHs 

In all the formulas to follow, the radius (a) and the wavelength (A) are in the units of 
10~^ cm unless stated otherwise. The choice of units is made to confirm with our earlier 
work (RSG-1 and RSG-2) wherein, the same units are used in the formulas for extinction 
spectra of the sihcate and the graphite particles. The advantage of using these units is to 
have the numerics associated with u {— 1/A) and a as coefficients or in radicals appearing 
in the various formulas are neither too large, nor too small. 

2.1.1 Absorption efficiency 

For absorption efficiency, we arrive at the following formulas: 
lOOOA < A < 1800A (FUV) 

Qabsi^) = X 

6.0 



0.93265 19.608, ^^,^,0 22.2817, „ 

+ (u - 0.643f (u - 0.7023f + 

u 1/ u 



^ ,i/- 0.7023)(1.0-i/)(i/-0.55) 



, (3.5A <a< 50A) 



where x — 2T:av. Note that the absorption efficiency is proportional to x. This is as it should 
be. For particles small compared to the wavelength of the incident radiation (x <^ 1), the 
general expression for the absorption efficiency, as a power series in is known to be of the 
form (see, for example, Pendorf [15]): 



Qabs{x) ^ x{A + Bx^ + Cx^ + ). 



(2) 



where A, B, C etc. are functions of the refractive index of the scatterer relative to the 
refractive index of the surrounding medium. The minimum and maximum values of x in the 
present study are 0.000977 (corresponding to a = 3.5 A ; A = 22,500 A ) and 0.314 (corre- 
sponding to a = 50 A ; A = 1000 A ) respectively. Therefore, in obtaining the formula (1) 
the first term of the general expression (2) is found to be sufficient in parametrization of the 
exact absorption efficiency successfully. Higher order terms in x give negligible contribution. 

Li and Draine [9] have given two formulas for the absorption efficiency to cover this spectral 
region which are as follows. For lOOOA < A < 1300A, 



Q 



abs 



100.0(1.286)% 
while, for 1300A < A < 1700A, 
a 



66.302 - 24.367i/ + 2.950i/^ - 0.1075i/^ 



(3) 



Q. 



abs 



100.0(1.286)% 



S2{X) + 1.8687 + 0.1905i/ + 0.4175(z/ - 5.9)^ + 0.04370(z/ - 5.9)^^ 



, (4) 



where 5*2 is the particular case {j = 2) of the general expression for Drude's spectral absorp- 
tion profile: 

i 'yjXjO'intJ 



7r(A/A,-A,/A)2 + 7: 



2 ■ 



(5) 



Values of 7j, Xj and aint,j are given in [9]. The expressions obtained therein are for absorp- 
tion cross section per C atom. The corresponding expression for the absorption efficiency is 
given in the appendix. 



1800A < A < 5000A (UV) 

In this wavelength region too, the parametrization of the absorption efficiency could be 
achieved in a form such that it is proportional to the size parameter x. 
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1.0 



0.74i/ + 83i/(1.0- ^)2_ 
The formulas obtained in [9] for this region are: 



, (3.5A <a< 50A) (6) 



Qahs — 



100.0(1.286)% 



S2{\) + 1.8687 + 0.19051/ 



(7) 



for 1700A < A < 3030A, and 

n 

Qabs — 



100.0(1.286)% 
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34.58 X 10-3-431/'^ cutoff(A, Ac) + 

i=3 



(8) 



for A > 3030A , where Sj is given by equation (5). It may be noted that cutoff (A, Ac) is a 
function of a as well as A because Ac is a function of the radius a of the PAH particle [9] . 

4000A < A < 8000A (Visible) 



(28.176i/^ - 0.091i/ - 0.311) - (a - 0.048)^ (4600i/^ - 2052.66i/ + 226.26) 
(0.041/100a)5/2 



l + (2175a)3 [0.1228 +M 



0.01312 _ 

2 



Figure 1(a) shows a comparison of absorption efficiencies computed on the basis of our formu- 
las with the exact results [14] for four PAHs sizes. The wavelength range is lOOOA to 8OOOA. 
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Figl (a): Comparison of predictions (lines) of absorption efficiencies from our equations with exact results 
(Draine website [8]) points 

(b): Comparison of predictions (lines) of absorption efficiencies form Li and Draine 2001 equations with 
exact results (Draine website[8]) points 

Frequency range 0.125 < < 1.0, wavelength range lOOOA < A < 8000A(i.e. FUV+UV+ Visible region) 




It can be seen that the agreement between the predictions of (1) and (6) and the ex- 
act results is excellent. To contrast results from formulas (1) and (6) with those from the 



parametrization by Li and Draine [9], a typical comparison of absorption efficiencies com- 
puted from formulas (3), (4), (7) and (8) with exact results is shown in Figure 1(b). Clearly, 
the accuracy of Li and Draine parametrizations and our formulas is almost identical. 

8000A < A < 12,500A (Infrared- I) 
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32.3a)" 


(0.00422/100a)2 
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- 1.0)2 
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(lOA <a< 50A) (10a) 



0.00208a 



1.0-2.4 



,+ f^l (.-0.08)(0.1245-.), 

0.22 + a(94.0 + C)](l.O-^) V « / 



(SA < a < lOA) 



(106) 



where 



C = 438(1.0 - 100a)(100a - 0.5)^. 



and the superscript s in (10b) stands for small PAH range SA < a < lOA. Figure 2(a) 
exhibits a comparison of variation of absorption efficiency computed on the basis of (10a) 
with exact results. For clarity of display, the comparison has been shown only for two values 
of PAH sizes, a — 20A and a — 50A. The agreement can be seen to be excellent. Although 
not shown here, computations have been done for (10b) also and the results have been found 
to agree reasonably well with exact results. 

It appears that the molecular structure of PAHs become vital in this wavelength range, for 
the absorption does not seem to follow the general expression (6) which is obtained from the 
expansion of Mie formulas for small x values. 



12,500A < A < 22,500A (Infrared- II) 

(SA < a < 50A) 
(11) 

A comparison of absorption efficiencies computed from (11) with exact results is depicted in 
Figure 2(b). Once again, the agreement of our formula with exact results is excellent. 
Above comparisons show that our formulas for absorption efficiency give excellent results 
over the entire wavelength and size range considered here. Formulas obtained involve only 
frequency of radiation and the size of the scatterer and are simple to use. In contrast. 
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Fig 2: Comparison of predictions (lines) of absorption efficiencies of our equations with exact result (Draine's 

website [8]) points 

(a) : Frequency range 0.08 <v < 0.125; wavelength range 8000A < A < 12,500A (IR-1) 

(b) : Frequency range 0.044 <v< 0.08; wavelength range 12,500A < A < 22,500A (IR-2) 

the Li and Draine formulas [9], for A > 3000A , require contributions from as many as 12 
Drude profiles making it somewhat cumbersome to use fluently. The formula also involves 
a function cutoff(A, Ac) containing the cutoff wavelength Ac which in turn, depends on the 
PAH radius in a complex manner. 

2.1.2 Scattering efficiencies 

For the scattering efficiencies, the formulas for the neutral or the ionic PAHs turn out to be 
identical. These are as follows. 



lOOOA < A < 1800A (FUV) 

QscaiX) — *^ 
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Fig 3: Comparisons of predictions (lines) of scattering efficiencies from our equations with the exact results 

(Drainc's website [8]) points 

(a) : Frequency range 0.55 < ly < 1.0; wavelength range lOOOA < A < 1800A (FUV) 

(b) : Frequency range 0.25 <u< 0.55; wavelength range ISOOA < A < 4000A (UV) 

(c) : Frequency range 0.044 <u< 0.25; wavelength range 4000A < A < 22,500A (Visible+IR) 

4000A < A < 22,500A (Visible, IR-I, IR-II) 

<3sca(3) = 1.66 x^ (SA < a < 50A) (14) 

The functional forms of the scattering efficiency in these spectral regions may be con- 
trasted with the expansion of exact Qsca in powers of x for small x (see, for example, Pendorf 
[15]): 

Qscaix) ^ x\D + Ex"" + Fx^ + ) (15) 

where D, E, F etc. are effectively, functions of the refractive index. It may be noted that 
while for (13) and (14) only the first term (of the order of x^) is sufficient to reproduce the 
observed results, for (12) the second and the third terms, of the order of x^, x"^ are also 
required. 



2.2 Ionic PAHs 



In the wavelength regions 1 and 2 (FUV, UV) the formulas remain the same as that for 
non- ionic PAHs. For the other three regions (Visible, IR-I, IR-II), following formulas have 
been obtained for the absorption efficiencies, 

4000A < A < 8000A (Visible) 



28.176i/^ - Omiu - 0.311 - (a - 0.048)^(4600i/^ - 2052.66i/ + 226.26)- 



(0.4858/lOOa)^ 



1.0 + 20.0(1001/ - 14.3)^ 
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1.0 + 0.20(100i/ - 13.4)' 



-, (lOA <a< 50A) 



(16a) 



0.0005(0.005/a) 



1.0 + 7(1001/2-1.96)2, 
8OOOA < A < 12,500A (Infrared- I) 

0.1a - 0.0001 



Qiabs (4) 



1.0 + [68.0 + ^ + 37.55(^ - l.o) ] - l.o) 



(5A < a < lOA) (166) 



r, (lOA < a < 50A) 
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Qtabsi^) 



0.21a - 0.001077 - O.OlflOOa - 0.768) 



+ 
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12,500A < A < 22,500A (Infrared- II) 

Qiabs{^) = Qa6.(5) + a [(0.005)2 ^ Q_;^^g(^ _ 0.0675)2) 



0.00352a - (0.04134 - a) 



022 ^089 _ ]^Q~j 



1.0 + 

(lOA < a < 50A) (18a) 

g|^^^(5) = a(0.3i/2 - 0.00434i/ + 0.000532), (5A < a < lOA) (186) 

where the subscript iahs used in Q indicates absorption efficiency for ionic PAHs. In the 
far-ultraviolet and ultraviolet regions, the absorption and scattering efficiencies for neutral 
and ionic PAHs are identical. They, however, differ in visible and infrared regions as can be 
seen from Figure 3 which depicts this difference for a = lOA and a = 50A . The curves have 
been plotted using above formulas. The curves show that the absorption by the ionic PAHs 
is generally greater in comparison to the neutral PAHs. More recently, Cecchi-PesteUini [16] 
have reported on variation of spectral properties of PAHs in various charged states. However, 



in the present work we limit ourselves only to the gross classification of neutral and charged 
PAHs, which is sufficient for our purposes. Scattering of electromagnetic radiation by a 
charged sphere has been treated by Klacka and Kocifaj [17] employing a generalization of 
the Mie theory. 

3 Extinction spectral features of PAHs 

Having obtained the absorption and scattering efficiency formulas in the various frequency 
ranges, one can generate the extinction spectra for PAHs for a specified particle concentration 
NpAH and a size distribution function fp{a) over a radius range [ao,am]- 

5A < ao, ttm < 

50A using the expression: 



^eM) = fj QLn^ym da, (19) 



fp{a) being normalized to unity. In (19), Q^^^{x) = Qabl^ + Q^^^ is the extinction effi- 
ciency of an individual scatterer of size parameter x — 271 a/ A; a and A being reckoned in unit 
10~^ cm. This choice of units is in keeping with our earlier formulas for extinction spectra 
of sihcate and graphite grains expressed in [1,2]. It may be observed that Qf^^ <^ Qabl^ 
the wavelength and radius ranges considered here. Hence, in (19) one can simply use Qab^^ 
in the place of Q^^^ without incurring appreciable error. 

The K^xt spectra (19) assumes much simpler form in the wavelength region where Qabt^ 
is linear in a. We observe that such is the case with the FUV and UV regions. In each of 
these regions we have Q^^^ ~ CL(f){i'). As a result, the corresponding expression for K^^f^ 
in the wavelength region lOOOA < A <5000A is given by the simple form 



- ^^PA^ <Pi-), (20) 



The volume corresponding to Npah number of particles will be 



VPAH = Y NpAH a^PAH X 10-l^ (21) 



where 



a^PAH = / a fp{a) da, 
a, ao, a^n all are in 10~^cm. units. Using (21) in (20), we finally have, 

K^rw « vpAH m- (22) 

In the other wavelength regions covering 5000A < A <22,500A , however, the K^^^ spectra 
would be dependent on the size distributional details in more involved manner rather than 



simply being proportional to the third moment as seen in (20). 

If a population of PAHs specified by Npah, fp{o) and Vpah is included in the carbona- 
ceous matter sector alongside graphite parallel as well as perpendicular components then 
due to stronger absorption features in the FUV and UV regions, a smaller amount of PAH 
(as compared to graphite) has the ability to increase the extinction contribution of the car- 
bonaceous component of the dust by much larger amount. As a result, when one has to 
include a population of PAHs to any sihcate-graphite model (two-component MRN model) 
in a manner so that the relative mass abundance criterion ( silicate mass/ carbonaceous 
mass) is maintained while the desired changes in the dust (PAHs -|- silicate -|- graphite) 
extinction spectra are obtained, one can concentrate mainly to the FUV and UV regions. 
Thus among many, a simple option while extending a two-component silicate-graphite model 
would be to replace a definite amount of existing graphite (by mass) with equal amount (by 
mass) of PAHs. Taking the mass density 2.24 gm/cm^) to be the same for both graphite 
and PAHs in bulk, we are thus led to a comparative study of the extinction spectra produced 
by equal volumes of graphite grains (parallel, perpendicular) and the PAHs. In the following, 
we examine this point in some details with respect to the MRN model considered by us in 



It may be recalled that in the two component graphite-silicate model of Mathis et al. 
(1979) [18], the spherical grains (graphite as well as silicate) are taken to follow a power-law 
size distribution of the form: 



where, a is the radius of the grain varying within the chosen values of ao and am- Besides, 
graphite material is taken to be present in two distinct structural varieties within the specified 
range of Oq and a^. This plausibility lies in the fact that graphite is a highly anisotropic 
material. The refractive index of graphite, therefore, depends on the orientation of electric 
field relative to its structural symmetry. Owing to practical difficulties in calculations of 
exact scattering (extinction) quantities due to the anisotropy, researchers have taken resort 
to an approximation known as "| — |" approximation [19]. In this approximation, graphite 
grains are represented as a mixed population of isotropic spheres, of which f» | fraction have 
refractive index m = m\\ (referred to as graphite parallel) and ~ | fraction have the refractive 
index m = m± (referred to as graphite perpendicular). By virtue of [1,2], we have formulas 
for generating graphite extinction spectra for both the components ( graphite perpendicular 
and graphite parallel) within the parameter ranges 



[1,2]. 




Gq < a < a 



(23) 



0.002//m < Oo < 005/xm; 0.15//m < < 0.25//m. 



By choosing Ngra — 4.4 x 10^ as the particle concentration for both varieties of graphite 
(perpendicular as well as parallel) and the size range [0.005//m : 0.25/xm], the corresponding 
extinction spectra are displayed in Figure 4(b). 
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Fig 4 (a): Comparison of predictions of absorption efficiencies of Neutral and Ionic PAH's. 

(b): Extincition spectra for equal volumes of PAH's and, parallel and perpendicular graphite grains. 



The associated material volume of each of the graphite components is Vgra = (47r/3) Ng^a a^graX 
lO^^^cc, a^^„ being the mean cube radius of the corresponding population of graphite parti- 
cles. The radius a has units in 10~^cm; a^gra is to be evaluated using the formula [1] 



gra 



5aQ(n 



with ttQ — 0.05; = 2.5 and — (am/cio) 



50. 



(24) 



If we have to generate the PAHs extinction spectra corresponding to an equal bulk of 
graphite (perpendicular or parallel) having this fixed volume Vgra we must use this value 
for VpAH occurring in (22). The resultant PAHs extinction spectrum is also displayed in 
Figure 4(b) along with that of the two graphite components bringing out the contrast in the 
individual extinction spectra produced by equal volumes of PAHs and the two varieties of 
graphite. It is clear from Figure 4(b) that the contribution of some mass of PAHs to the 
interstellar dust extinction will be much greater in comparison to equal mass of graphite in 
the wavelength range lOOOA < A < 5000A (FUV and UV regions). Further, if one looks at 



the ratio Kf.xt{^ / Kext{\) fixing at 5000A and compares among the three species viz. 
PAHs, graphite (perpendicular) and graphite (parallel), one can see easily that the ratio is 
much greater for PAHs as compared to that of the other two species within the wavelength 
range lOOOA to 5000A. Clearly, replacement of a small amount of graphite by an equal 
amount of PAHs in a dust mixture will cause an increase in the dust extinction throughout 
the wavelength range lOOOA to 5000A. However, this increase would be greater in the FUV 
region than that in the UV region. This is how the addition of PAHs to sihcate-graphite 
model will help one to adjust the interstellar dust extinction spectral features in the FUV, 
UV regions. 

In connection with our choice for the PAHs grain size range, we would like to mention 
that, in Li and Draine [9] size range of PAHs has been taken to be 0.00035 < a < 0.01 jim 
in their investigations. The absorption cross section of carbonaceous grains is represented 
as [9] 

CZ\a, A) = ePA^^Cr,^^(a, A) + (1 - iPAH)Ci:{a, A), 

^PAH{a) = (1 - qgra) X mm[l, {a^/af], 

= 50A , qgra = 0.01, 

where C^^^ ^ C^l^ are, respectively, absorption cross-sections of PAH and graphite grains 
of radius a and wavelength A; is the grain radius from which the transition from PAH 
properties to graphite properties begins; ^pah is the PAH weight factor which falls off con- 
tinually as a increases from to infinity. However, as has been mentioned in [9], this 
choice of the weight factor is quite arbitrary. With this freedom of choice for C,pah, we con- 
sider a dust model in which the carbonaceous component is made up of PAHs of size range 
5A < a < 50A and graphite grains in the size range 0.005/xm < a < 0.25/xm. 

4 Conclusions and discussions 

In this work, we have analyzed the extinction spectra generated by small size carbonaceous 
matter present in interstellar dust in the form of Polycyclic Aromatic Hydrocarbons (PAHs). 
Our investigations cover the wavelength region lOOOA < A <22,500A with PAHs size range 
SA < a <50A , a being the radius of the PAH particle. In doing so, we have made use of 
the scattering and absorption efficiency data available in [14] to obtain simple and accurate 
analytic formulas for scattering and absorption efficiencies expressible as functions of a and 
A. These formulas can be used expediently for computing extinction spectra generated by 
some bulk PAHs material specified through any selected size distribution. Thus, this work, 
together with our earlier work [1,2] essentially completes our search for an analytic platform 
for interstellar dust extinction spectra analysis which employ MRN type dust models con- 
taining sihcate, graphite and PAHs as well. 



The scattering efficiency has {a/Xf behavior over I8OOA < A <22,500A. For lOOOA < 
A <1800A, the behavior is shghtly different in that additional {a/Xf and {a/Xf behavior 
also creeps in. This is quite in keeping with the Pendorf formula [15] for the scattering 
efficiency in the case of small values of a/ A. The absorption efficiency has a/ A behavior 
over lOOOA < A <5000Awhich is again in keeping with the Pendorf formula for absorption 
efficiency in respect of small values of a/A. However, for A >5000A sharper absorption fea- 
tures characterizing the special dispersive properties depending on size of the particle appear. 
This results in the absorption efficiency having a functional form which differs from that of 
Pendorf. A good account of modelling for these sharper absorption features occurring in the 
extinction spectra data experimentally generated by small sized hydrocarbons prepared in 
the laboratory is given by Li and Draine [9] . 

The ionic PAHs show a shghtly different extinction behavior as compared to their neu- 
tral counterparts. Both the neutral and the ionic species have the same scattering efficien- 
cies and hence scattering properties of these are exactly the same within the size range 
3.5A < a <50A and the wavelength range lOOOA < A <22,500A. The neutral and the 
ionic PAHs absorption efficiencies are same for lOOOA < A <4000A; but they differ in the 
wavelength range 4000A < A <22,500A. The pattern of variation is illustrated in Figure 
4(a). One can see that the range of wavelength [Ai, A2] over which the absorption efficiency 
differs depends on the PAH radius a. It appears that both Ai and A2 increase with a. As 
the efficiency differences are not substantial, the preference of ionic over neutral PAHs to be 
included in the dust model does not seem to provide with much changes in the corresponding 
extinction spectra. 

To estimate the extent to which the extinction spectra of a dust can be affected by hav- 
ing a definite amount of PAHs in it, we have compared the extinction spectra generated 
by equal volumes of all the three carbonaceous components viz, graphite (perpendicular), 
graphite (parallel) and PAHs. As PAHs extinction is essentially due to absorption only, 
we have chosen the wavelength range 1000 A < A <5000A due to the simple form of the 
absorption efficiency. The comparison brings out very clearly the fact that PAHs indeed can 
act as a very effective carbonaceous component in MRN type dust models. It will be our 
future endeavor to use the analytic framework which includes this work alongside [1,2] for 
the analysis of extinction spectra data corresponding to MW, LMC and SMC. 



Appendix 



The absorption cross section per C atom has been denoted in [9] as C^i^^ /Nc, where Nc 
is the number of C atoms. The absorption efficiency, Qa^^i for a PAH of size a can be 
expressed as 

The "radius" a of a PAH containing Nc number of C atoms is defined to be the radius 
of a sphere with carbon density of graphite containing the same number of C atoms, i.e., 
a = 1.286A^c^l Thus, 

where C^f^/Nc is given in equations (5)- (11) in [9] for various spectral ranges. 
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